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^ , Abstract 

We investigate the quantum mechanical wave equations for free 
(3JT), particles of spin 0, 1/2, 1 in the background of an arbitrary static grav- 

itational field in order to explicitly determine if the phase of the wave- 
function is S/Ti — f p^dx^/h, as is often quoted in the literature. 
' We work in isotropic coordinates where the wave equations have a 

simple managable form and do not make a weak gravitational field 
approximation. We interpret these wave equations in terms of a quan- 
tum mechanical particle moving in medium with a spatially varying 
effective index of refraction. Due to the first order spatial derivative 
structure of the Dirac equation in curved spacetime, only the spin 1/2 
particle has exactly the quantum mechanical phase as indicated above. 
The second order spatial derivative structure of the spin and spin 1 
wave equations yield the above phase only to lowest order in Ti. We 
develop a WKB approximation for the solution of the spin and spin 
1 wave equations and explore amplitude and phase corrections beyond 
the lowest order in %. For the spin 1/2 particle we calculate the phase 
appropriate for neutrino flavor oscillations. 



'Contact author for correspondences. 
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1 Introduction 



The phase of a quantum mechanical particle in curved spacetime has been of 
considerable interest both theoretically and experimentally for many years. 
In the late 1970's researchers were interested in explaining the quantum 
mechanical interference fringes for neutrons traversing different paths in the 
Earth's gravitational field 0]. Recently, there has been renewed interest in 
this topic in order to investigate the interplay between gravitation and the 
quantum mechanical principle of linear superposition in relation to flavor 
oscillations of neutrinos in the context of the solar neutrino anomaly, and 
type-II supernova §2|, ||, §[]. 

Most computations of the quantum mechanical phase (QMP) of a free 
particle in curved spacetime refer back to the seminal article by Stodolsky 
H who argued that the relativistically invariant phase S/% would be given 
by 

S(r,t)/h = - / ppdaf. (1) 

n Jr A ,t A 

In Eq. (jl|) = mg liu u 1 ' is the general relativistic 4-momentum, and = 
dx^/dr is the 4- velocity such that g fiU u ll u L ' = 1. The phase is computed 
along a geodesic connecting the points [ta, £a) and (r#, ig). The quantum 
mechanical particle is assumed to be a test particle in the sense that it moves 
in the background of the general relativistic metric and does not generate it 
own gravitational field. For any arbitrary metric in general relativity (GR), 
static or otherwise, we have a relationship amongst the momenta, called the 
mass shell constraint, given by 

d^VvVv = (mc ) 2 , (2) 

where m is the rest mass of the particle and cq is the vacuum speed of 
light. The form of the scalar quantum mechanical wave function proposed 
by Stodolsky was 

i/)(r,t) =Ae iS{r ' t)/h , (3) 

where the amplitude A is assumed constant. 

It is clear from Eq. (|3|) that a single spatial derivative of ip(r, t) can 
generate p^, but that two spatial derivatives cannot generate the mass shell 
constraint Eq. (||), since imaginary cross terms involving gradients of p^ 
will be produced. Therefore, for second order wave equations, for which 
spin and spin 1 particles are particular cases, the QMP cannot exactly 
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take the form of Eq. (|l]). In addition, the form of ip(r,t) in Eq. (||) must 
involve an amplitude change, which is well known from the standard WKB 
approximation of the wave function in non-relativistic quantum mechanics 

In this paper we examine the general relativistic wave equations for quan- 
tum mechanical particles of spin 0,1/2,1 in the background of a generic 
static metric of arbitrary strength which can be written in isotropic form. 
We do not make a weak gravitational field approximation, which is often 
the case in the literature By working in isotropic coordinates, the wave 
equations are expressed in simple forms, which are then easily intrepreted. 
We explicitly demonstrate that for the Dirac equation in curved spacetime 
the form of the QMP is given exactly by Eq. (jl]) , and that Eq. (Q) is exactly 
satisfied. This result is directly traceable to the first order spatial derivative 
structure of the Dirac equation. The only difference from Stodolsky's pro- 
posal is that amplitude in Eq. (||) takes on a simple spatially varying form. 
For particles of spin and spin 1, the wave equation is of second order in 
both time and space. We develop a WKB approximation for ijj(r, t) which to 
lowest order in h has the form of Eq. (Q) and satisfies Eq. (^). We formally 
solve for if) to all orders in fi. However, since the WKB expansion is only an 
asymptotic series § we only examine the next higher order phase and am- 
plitude corrections. We give an interpretation of the above wave equations 
in terms of a particle moving in a medium with a spatially varying index of 
refraction n(r) as discussed in Evans et. al. || and Alsing [[H]] (see also 

0)- 

The outline of the paper is as follows. In Section || we briefly review 
Stodolsky's reasoning for the form of the phase of a quantum mechanical 
particle in curved spacetime as given in Eq. (||). In Section || we examine the 
general form of the quantum mechanical wave equation for a free particle in 
a background curved spacetime. In Section |I] we examine the Dirac equation 
in curved spacetime and show that in isotropic coordinates, it takes on a 
simple form even for arbitrary strength gravitational fields. In Section || 
we examine the wave equation for a scalar particle in curved spacetime and 
develop a WKB solution. We develop ^-dependent quantum corrections 
beyond the lowest order WKB amplitude and phase approximations found 
in standard textbooks. This WKB approximation has applicability to the 
curved spacetime spin 1 wave equation since the latter has the same form 
as the scalar wave equation for each of the vector components of the wave 
function. Throughout the discussion we draw the analogy to a particle 
moving in a medium with a spatially varying index of refraction, which the 
form of these wave equations explicitly exhibit. In Section |6| we relate the 
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classical optical-mechanical analogy of the lowest order approximation to the 
QMP. In Section g we calculate the QMP for a spin 1/2 particle appropriate 
for neutrino flavor oscillations in a gravitational field and relate it to the 
effective index of refraction discussed in the previous section. In Section || 
we summarize our results and discuss their implications in light of current 
calculations in the literature. 



2 Stodolsky's proposal for the phase of a quantum 
mechanical particle in curved spacetime 

In 1979, Stodolsky |J argued the phase of a spinless quantum mechanical 
particle should take the form of Eq. ([l]). In the following we summarize 
his reasoning. In flat Minkowski spacetime elementary quantum mechanics 
tells us the phase of the particle's wavefunction is given by dimensionless 
quantity (p-x — Et)/h. The numerator of this expression has units of action, 
mcox (proper distance). Since the particle is taken to be a test particle, 
and therefore does not generate its own gravitational field, classically it 
would follow a geodesic. This is assumed to be true quantum mechanically 
as well, independent of the particle's spin. For an arbitrary metric ds 2 = 
g^y dx^ 1 dx v , Stodolsky proposed the phase <& accumulated by the particle in 
traversing the geodesic connecting the spacetime points A and B to be 

$ = ^ [ mds = S/h. (4) 
n J a 

The particle's quantum mechanical wavefunction is then taken to be pro- 
portional to the phase factor e l *. The integral appearing in Eq. (||) is just 
the relativistic action for a particle moving on a geodesic JE|. Stodolsky 
states (||, p392), "If we wish the phase to be an invariant and to agree 
with elementary quantum mechanics this seems to be the only choice." Fur- 
thermore, dividing the metric by ds and defining = thcq g ^ dx v j 'ds as 
the canonical general relativistic momentum yields Eq. (II). Stodolsky finds 
additional reassurance from the observation that Eq. (|l) also appears rea- 
sonable from the Feynman path integral approach to quantum mechanics. 
If Eq. (||) is valid for paths neighboring the classical path then one obtains 
in the usual way that the actual classical path is the one for which (5$ = 0. 
But from the form of Eq. (||), this is just the classical relativistic condition 
for a geodesic. In his paper Q Stodolsky examines Eq. (||) in the limit of 
weak gravitational fields for a for both static and stationary metrics. 
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The purpose of this work is to explicitly test whether Eq. ([l]) is a valid so- 
lution of the general relativistic wave equations for particles of spin 0, 1/2, 1. 
To this effect, we consider an arbitrary static gravitational field which can 
be written in isotropic coordinates Jl^j . For such gravitational fields, the 
spatial portion of the metric can be written in a conformally flat form, 

ds 2 = n 2 (r) cl dt 2 - (dx 2 + dy 2 + dz 2 ) = n 2 (r) cl dt 2 - — 1~- \dr\ 2 . 

(5) 

In the above, r is the isotropic radial marker. 

Many static gravitational fields of physical interest can be written in the 
form of Eq. (||). The important class of static, spherically symmetric metrics 
written in the usual spherical coordinates {t, r' , 9, <p} have the general, non- 
isotropic form 

ds 2 = ^r ) c odt 2 -^r ) - R\r') ( de 2 + » 2 (9) # 2 ). (6) 

For example, the Schwarzschild metric has F{r') = (1 — r' s /r'), G(r') = 1, 
and R(r') = r' where r' is the usual radial coordinate measured by an 
observer at infinity and r' s = 2GM/cq is the Schwarzschild radius. After 
transforming the Schwarzschild metric to isotropic coordinates {t, r, 9, 4>} one 
obtains ^(r) = (l+r s /r)/(l— r s /r) and <3?(r) = (l+r s /r) -2 where r s = r' s /2. 
The relationship between the coordinate radius r' and the isotropic radius r 
is given by r' = r<I>~ 1 (r) = r (l + r s /r) 2 . Other interesting static spherically 
symmetric metrics for both Einstein gravitation theory and non-Einstein 
gravity theories can be found in |T4| . 

The utility of isotropic coordinates stems from the ease by which we can 
define an effective index of refraction for the gravitational field. By setting 
ds = in Eq. (||) we obtain the coordinate speed of light c(r) as 



c(r) 



dr 



dt 



c n(r)<Z>(r)^^- (7a) 



n{r) s TMWY (7b) 

The paths of both massive and massless particles in such metrics (i.e. geodesies) 
can be interpreted as motion through a medium of index with an effective re- 
fraction n(r). For the Schwarzschild metric we have n(r) = (1 + r s /r) 3 /(l — 
r s /r). Evans et al. M used this concept to write the geodesic equations 



of motion for static metrics in Newtonian "F=ma" form. Alsing [1C] later 
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extended this idea to the case of stationary metrics. In this paper, the quan- 
tity n(r) will continually arise in the wave equations developed, and will be 
naturally interpreted as an effective index of refraction. 

For the metric given by Eq. (||) we can write the mass shell constraint 
Eq. © as 

where we define the 4- momentum as = (E/co, p). This expression re- 
duces to the ordinary flat space, special relativistic form E 2 = (pco) 2 + 
(?h,Cq) 2 as r — > co. This naturally leads us to propose the quantization rules 

E = huj, p = hk, (9) 

which define the frequency & and wave vector k. Using Eq. (||) in Eq. (g), 
we can rewrite the mass shell constraint as 

n 2 (r)k 2 +k 2 (r) = k 2 c (r). (10) 

In Eq. ( JT0| ) we have defined the 4- wave vector as k^ = (kg, k) with kg = uj/cq, 
k 2 = k k and k c {r) = {\ c <&{r))~ l where A c = %/ (m Co) is the usual Compton 
wavelength for a particle of mass m. Note the wave vector k(r) and hence 
the velocity of the particle is position dependent as appropriate for a particle 
in a medium with index of refraction n(r). 



3 Quantum wave equations in curved spacetime 

In practice, we only know how to quantize wave equations in flat Minkowski 



space 1 15]. Therefore, a natural way to describe curved spacetime is to 



erect local coordinate axes named vierbeins or tetrads [16, 17, 18] at each 
point X in spacetime and then project all tensor quantities onto these local, 
Lorentzian inertial frame axes. At each point X the local metric takes the 
flat spacetime form, r] a f, = diagonal {1,-1,-1,-1}. Throughout this paper, 
Latin indices {a, b, c, . . . } near the beginning of the alphabet will refer to the 
local inertial frame with values {0, 1, 2, 3} while Greek indices {^u, v, A, . . . } 
will refer to the general coordinate system with values {0, 1,2,3}. The 
tetrads e fl a (x) and e a M (x) are defined by 

9^(x) = e-{x)e v \x)^ ab , gT{x) = e a »(x) e b »{x) V ab . (11) 



6 



Contravariant and covariant vectors and in the general coordinate 
system can be expressed as vectors V a and V a in the local inertial frame 
(and visa versa) by means of the transformations 

V»{x) = V a (x) e/(x), V a (x) = V»(x) e/(x), 

V^x) = e°(x) V a (x), V a (x) = e a »{x) V^x). 

Local inertial frame indices are raised and lowered with the flat spacetime 
metric r] ab , while all general coordinate frame indices are raised and lowered 
with the metric g^. 

The prescription for generalizing a flat spacetime wave equation to curved 
spacetime proceeds as follows: (i) begin with the appropriate flat spacetime 
Lagrangian for the wave equation of interest, (ii) replace all local partial 
coordinate derivatives by covariant derivatives via d a — > e a ^(x) V^, and (iii) 
contract all vectors, tensors, etc. into n-biens, (V a (x) — * V^(x) e fl a (x) , 
etc.). 

For a field ip(x) of arbitrary spin the spin covariant derivative V M is 



defined by |§ [TJ, U> M 



V u ifj(x) = [d v + tt v (x)]^{x), (12a) 
Q u (x) = - % -u abu (x)o ab = \uj abv (x)[l a ,l b l spin 1/2 (12b) 

£lv{x) = 2 ^abuix) T, ab . arbitrary spin (12c) 



In Eq. ( |12|) the Fock-Ivanenko coefficients Q. v {x) (not to be confused with 
^( r ) — v 9oo( r ) m ^ ne isotropic metric Eq. (|5|)) are defined in terms of the 
spin connection coefficients 0J a bv given by 

= S° [%%u = e f (due^+e^T^), (13a) 
^abu = e^e^, (e^ a = g^ u rj ab ej^j , (13b) 

which are antisymmetric in their Lorentzian indices 0J b av = ~~ ^abv The 
semicolon denotes the usual Riemannian covariant derivative, V^- v = d v V^ + 
T^wV 7 . Since r] ab;a = and the Riemannian metric compatibility condition 
gives g^v-o = 0, we can freely raise and lower both Lorentzian (Latin) and 
general coordinate indices (Greek) within a Riemannian covariant derivative 
operation. Note Eq. (|l3a| ) can be rearranged to define the action of the spin 
covariant derivative V u on the tetrad e a M , 

V„ e a " = d u e/ + e a CT - J av e b » = 0, (14) 
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analogous to the Riemannian metric compatibility condition. Mnemonically, 
for the spin covariant derivative of a quantity T^'"^'" with a mixed set of 
index types, each general coordinate index receives a contribution from the 
metric connection and each local Lorentzian index receives a contribution 
from the spin connection. 

Finally, in the last expression Eq. ( p. 2d ) the constant matrices T> ab are 
the generators of the Lorentz group for an arbitrary value of the spin |l6| , 
18 1 . For the specific case of spin 1/2, S afe = l/4[7 a , r y b ], where r ) a are 



the constant Lorentzian gamma matrices, with an explicit representation 
given in Appendix [A]. In addition, we have used the conventional notation 
a ab = i/2[j a , 7 fe ] in Eq. (ggg ). 

We are now ready to generalize the flat spacetime Lagrangian C(x) for 
the fields of interest to curved spacetime. Since in flat spacetime the action is 
given by Sji at = J d 4 x£(x), the curved spacetime analogue of C{x) must in- 



volve the volume factor y/—g{x) = det^e M a (x)J , where g{x) = detyg^(x] 
in order that S curve d transforms as a scalar. 

3.1 Spin field 

For a spin field <f>{x), T, ab = 0, d a <f){x) -» e/V M 0(z;) with V ^(x) 
dfj,4>(x), and we have 

2 „2 

o 



C(x) = \ V ab d a cf>d b 



2 \ ' ur UT h 



2 



- ~(-g) 1/2 (ri^d^e^d^-^A. (15) 

Variation of the Lagrangian with respect to <j)(x) yields the general relativis- 
tic Klein-Gordon equation (GRKGE) 

( □ + ^ ) m = <T <t>w + (mc /h) 2 <i>(x) 

= (^d^-g^I%d x + ^)tf>(x) = 0, (16) 
where we have defined the covariant Laplace-Beltrami operator □ 4>(x) = 



L. v . For the case of spin we can actually derive Eq. (16) much more 



directly from Eq. (|2|). Classically, the Lagrangian for geodesic motion is 
given by L = m cq ds = m cq w 'g^ u x^ x u , where x^ = dx^/dr, dr = ds/cQ. 
The Hamiltonian is given by £^ (dL/dx^) x^ — L = 0. This implies there 
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must be a constraint amongst the momenta, which is given by Eq. (0), 



and acts as the effective Hamiltonian for the system 21], H e ft = g^ v p^Pv 
Substituting in the quantization condition p^ — ► —ifrV^ yields Eq. ( |l6| ) 
directly. 

Note the most general Lagrangian allows for a term proportional to the 
Ricci scalar R(x) = R^(x), i.e. an extra term R(x) cj) 2 (x) in Eq. (|i~5|), 
which leads to a corresponding term £ R(x) <fi{x) in the GRKGE |{T(f| . The 
value of £ is not determined by any physical principle and to date has not 
been experimentally measured due to the minute effects of curvature in 
our solar system. The value of £ = is called minimal coupling for ob- 
vious reasons. In n dimensions the value of £(n) = 1/4 (n — 2)/(n — 1), 
which assumes the value of £ = 1/6 in four dimensions, is called confor- 
mal coupling. For massless particles under conformal coupling the GRKGE 
is invariant under the simultaneous conformal transformation of the metric 
9fii/(x) — > f 2 (x) g pLV {x) and field transformation 4>{x) — > f~ l (x) <j){x) for an 
arbitrary function f(x). In this paper we concern ourselves with regions of 
spacetime devoid of matter so that Einstein's equations reduce to R^ v = 
and consequently R{x) = 0. 

3.2 Spin 1/2 field 

For the case of spin 1/2 we have S ab = l/4[7 a , 7 fe ] = -i/2a ab , d a i\) 
ea^^fiip, and we obtain 

C(x) = l - (^7°^ " (0a$W " 

- (-9) 1 ' 2 ~ (jhV V M V - e/ (V^) 7 a V - (17a) 

= (-9) 1 ' 2 \ V M V - (V^) 7 ^V " ^</>) (17b) 

In the above we have defined the curved spacetime counterparts to the 
Dirac matrices as 7^(x) = j a e a ^(x) which satisfy the Clifford algebra 
{j^(x), 7 iy (x)} = 2g^ u (x). Variation of Eq. (|17b| ) with respect to ^(x) = 
ip^(x) 7°(x) yields the general relativistic Dirac equation (GRDE) for spin 
1/2 particles 

(i7^(x)V^-^)#E)=0. (18) 
Note that no geometric curvature term with a dimensionless coefficient can 



be added to the Lagrangian and thus to the GRDE [16] 
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3.3 Comment on the Pauli-Schrodinger equation 

In flat spacetime, the Dirac wavefunction also satisfies the Klein-Gordon 
equation due to the factorization (i~f a d a + mco/h) {i^ a d a — mco/h)ip = 
which implies \r] ab d a db + (mco/h) 2 ~\ip = 0. However, in the presence of 
a gravitational field, this is no longer the case |p2| . Using Eq. (|l8|), the 
analogous calculation produces (^7 M V M + mco/h) (ry^V^ — mco/fuip = 
=>• [VVV^V,, + (mc /h) 2 ]ip = 0. This can be written as [g^V^Vy - 
l/2a^K^ + (mco/n) 2 ] V = 0, where = 1/2(V„ V M - V„ V„) = - 
dpQ'v + [^i/)^u] is called the spin curvature in analogy to the Riemann 
curvature |23|]. The relation between the spin curvature and the Riemann 
curvature is given by (see appendix of |p2f ) K^ v = \/ <iR^ va pa a ^ and RI = 
Ri±vai3&^ v = —2K iW o ilv where R is the Ricci scalar and I is the 4 x 
4 unit matrix. Using the trace of Einstein's equations R^ u — l/2g^ u R = 
— 8-kG / CqT^ we obtain 



IttC 

g^V^V„ + — T+(mco/h) 2 
c o 



V> = 0, (19) 



where T = T^ is the trace of the energy-momentum tensor. 

Eq. (19) is the generally covariant extension of the Pauli-Schrodinger 
equation p4j and describes spin 1/2 particles in a gravitational field. As 



such, the covariant derivative in Eq. (19) is given by Eq. (12a) involving 
the Fock-Ivanenko coefficients. On the other hand, Eq. ( |i~6| ) is the generally 
covariant extension of the Klein-Gordon equation for spin particles, with 
no 0,^ terms. If we expand the spin covariant derivatives in Eq. ([[!]) using 
Eq. (|12a|) we have 

[•T^/A + (mc /h) 2 } iP + (2TtG/cl) 2 Tij 

+ g ^(2n lx d u + (d lx n u ) + n lx n u )^ = o. (20) 

In flat spacetime and in regions where T = 0, the tetrads are constant unit 
vectors and g^ u — > rf v in Cartesian coordinates. Hence, the Fock-Ivanenko 
coefficients are zero. The remaining terms in the square brackets above 
constitute the usual Minkowski Klein Gordon equation for the spinor field 



ip{x). Similarly, for the scalar wave equation Eq. (16) in flat spacetime in 
Cartesian coordinates, the affine connection T a ^ v is also zero, and so again 
we recover the Minkowski Klein Gordon equation, for the scalar field <j){x). 
The point to be made here is that in flat spacetime the Dirac wave function 
(which describes spin 1/2 particles) is also a solution of the Klein-Gordon 
equation (which describes spin particles) by the virtue of there being no 
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gravitational field p^ ]. In curved spacetime, this property no longer holds. 
In the presence of a gravitational field, the solution to the generally covariant 
Dirac equation is a solution of the generally covariant Pauli-Schrodinger 
equation. 



3.4 Massless spin 1 field 

In Minkowski spacetime the electromagnetic (massless, spin 1) field F ab 
is described by the Lagrangian C = -l/4F af, F a j. where F a t = d b A a — 
d a Ab in terms of the vector potential A a . The gauge freedom of the vector 
potential prevents a straightforward quantization of the theory and leads to 
the introduction of gauge-fixing terms Cq = —l/2(~ 1 (d a A a ) 2 where £ is a 
parameter determining the choice of the gauge. In the Feynman gauge £ = 1, 
variation of the combined action C + Cq with respect to the vector potential 
leads to the wave equations OA c = n ab d a d b A c = for each component of 
the vector potential. 

In curved spacetime the electromagnetic field takes on the same form 
as in Minkowski spacetime due to the cancellation of the connection terms, 
F^ v = A^ v — A UII1 = d u A^ — d u A u . To generalize the Maxwell field to 
curved spacetime we make the substitution A a — > A^ d a — > e^V^ with 
Qp(x) given by Eq. fll2c] ) with [T, ab ] c d = 5 a d r] bc - 5 b d rj ac corresponding to 
the (1/2,1/2) representation of the Lorentz group. The action takes the 
form of C = —l/AF^ F^ and correspondingly the gauge fixing term takes 
the form Cq = — 1/2 C _1 (V^^) 2 . Variation of the action C + Cq leads to 
the wave equations 

+ R»A v -(l-C X )AJ ^ = 0, 
□^4^, = 0, Feynman gauge and R^ u = (21) 

For our purposes, each component of the vector potential A^x) in Eq. ( |2l| ) 
satisfies the massless GRKGE, so that it suffices for us to examine the spin 
GRKGE Eq. ©. 



4 Plane Wave-like solutions for the GRDE 



We are now interested in finding solutions to the GRDE of the form Eq. (0) 
which we write as ip(r,t) = A{r ,t) e 1 s ^ v,t ^ h where u is a Dirac spinor. 
Writing out Eq. (|l8| ) and using Eq. ( |12b| ) the GRDE is explicitly given as 



-a ab e 



e-bv 



■41 



mco 



tp(x) = 0. 



(22) 
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Using the explicit representation of the Minkowski spacetime Dirac matrices 
7° given in the Appendix [A], substituting in the static metric Eq. (||) and 
dividing the resulting equation by $(r) yields after some algebra p3] 



i ( n(r) 7 °a + 7 • V + ^ 7 • V((x) ^ - k c (r) 



i>{x) = 0. (23) 



In Eq. (^) we defined V = {d/dx,d/dy,d/dz} as the spatial gradient in 
the general coordinate frame and 7 = {7 1 , 7 2 , 7 3 } the spatial components of 
the constant Dirac matrices in the local Lorentz frame. The quantity 
is given by C( x ) = hi[$7(r)/<]? 2 (r)]. 

VC(x) can be removed by defining 
V,t). A simple calculation reveals 
The resulting equation for 4>(x) is 



The middle term in Eq. ( p3|) | 7 
the wavefunction as ip(r,t) = f(r)4 
f(r) = e" 3 / 4 ^) = ($ 2 (r)/f)(r)) 3 / 4 . 
given by 



n(r) d 

7° — 757 + 7 • V 
cq at 



(r,t) = 0. 



(24) 



Eq. ( |24D is directly interpretable as the Dirac equation in flat spacetime with 
a spatially varying index of refraction n(r) and a spatially varying Compton 
wavelength A c (r) = \[ lat §(r), k c (r) = 1/A c (r), \[ lat = TijmcQ. 

We now seek a plane wave-like solution by substituting in Stodolsky's 
suggestion for the QMP in the form 



(r, t) = A exp 



^ ^ p(x) -dx - Et 



(25) 



where A is a constant Dirac spinor. Note that E is a constant of the geodesic 
motion given by E = po c o = i£l 2 (r)mc^ since the metric is independent 
of the coordinate t. Separating Eq. (|24] ) into a coupled set of two-spinor 
equations we have 



n(r) E/co ■ 
—a 



- mco/$(r) 
p(r) 



—a • p(r) 

n(r) E/cq + mco/&(r) 



A, 
A. 



0, (26) 



where a are the usual 2x2 Pauli spin matrices (Appendix [A|) . Eq. ( |26"| ) is a 
homogeneous set of algebraic equations and therefore only has a nontrivial 
solution if the determinant of the coefficients is zero. Using the identity 
(<r • p) 2 = p- p = p 2 and setting the determinant of the coefficients equal to 
zero yields precisely the mass shell constraint in the form of Eq. (||). 
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To summarize, we have found an exact plane wave-like solution for the 
GRDE for static metrics in isotropic coordinates to be of the form 



^(r,t) = A[ 1 ±±) exp 
d> 2 (r)^ /4 



~ ( J p{x) -dx - Et) 



= A [m) exp lv Kdx i- (27) 

This solution is normalizable with respect to the inner product 

rfV^d^^M^M), (28) 



where det gij = $ 6 (r) is the determinant of the spatial portion of the 
metric g^ u . 

Eq. (^?|) is one of the main results of this paper. It is an exact solution, 
valid for arbitrary strength gravitational fields for both massive and massless 
fields. Several authors g, [2§, |27|, f2j§ || begin with Eq. (|22|) but examine the 



equation in the weak field limit e a = S a + 1/2 h a . In some of those works, 



a Foldy-Wouthuysen transformation [[3C], 31] is performed to write the GRDE 
in the form of an effective weak field Schrodinger equation, ih dip/dt = H ip. 
Typically, one is interested in the Hamiltonian with respect to the measure 
J d 3 r (versus the measure in Eq. (p8|)) so that the momentum operator 
can be interpreted in the usual flat spacetime form as a spatial gradient in 
Cartesian coordinates, p = —ifi'V. In such cases, a new wave function is 
defined by ip = (det gij) 1 ^ 4 "^ = Uip, with the corresponding Hamiltonian 
H = UHU^ 1 . The focus of the above works is for the most part, the 
physical interpretation of terms in the effective weak field Hamiltonian H as 
post-Newtonian corrections to the gravitational potential, with an emphasis 
on spin-gravity coupling. 

Instead, in this work, we are interested in the solution of the GRDE 
Eq. ( p2|) directly, with emphasis on strong gravitational fields and the is- 
sue of whether or not the phase of the wavefunction takes the form pro- 
posed by Stodolsky. If one were to try force an interpretation of Eq. (|24|) 



as a Schrodinger equation, the effective Hamiltonian would contain non- 



Hermetian terms. This fact is well known (see [26, g7|) and arises because 
of our removal of terms that involved gradients of the gravitational field, 
and the introduction of a spatially varying function n(r) in front of the first 
order time derivative. As stated earlier, we interpret Eq. ( ^4|) as the flat 
spacetime Dirac equation for a particle moving in a medium with a spatially 
varying effective index of refraction and Compton wavelength. Regardless 
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of any imposed physical interpretation, the point we wish to emphasis here 
is that the GRDE admits a solution with exactly the phase suggested by 
Stodolsky Eq. (|I]), and a spatially varying amplitude. 



5 WKB Solution to the GRKGE 



To develop a solution for the GRKGE Eq. (Jlq) , we begin by using the 
identity @ fl^T*, = - {-g)~ 1/2 d^[{-g) 1 / 2 g^\. Substituting this and the 
static metric Eq. (B) into Eq. (]T6|) yields 



2 o2 / 

^^-V 2 V>-V£-VV> + A; 2 (r)V> = 0, (29) 

where we define the quantity £ = ln(J2(r)/$(r)j . As in the previous section, 
we can remove the term linear in by seeking a solution of the form 

V>(M) = /(r)0(r,i), (30) 

with f(r) = ($(r)/ri(r)) 1 / 2 and (j)(r,t) satisfying 

" W + [*2 ( r ) + ^ H <^ = °> ( 31 ) 



where 

??(r)^V 2 £(r) + i|V£(r)| 2 . (32) 

Note Eq. ( |3l| ) is essentially the flat spacetime wave equation for a scalar 
particle moving in spatially varying index of refraction and Compton wave- 
length. However, there is an additional term ij(r) in the wave equation 
which is a direct result of the geometric cross term g^ v T^ v arising from the 
product of the covariant derivatives. In the weak gravitation field limit this 
term is typically dropped |5| when seeking an approximate solution to the 
wave equation. This term can also be exactly eliminated by transforming to 
harmonic coordinates (|]l8| pl61-163) for which the identity above becomes 
four coordinate conditions, namely d^[{— fl , ) 1 ^ 2 fl r/iA ] = 0. However, the use 
harmonic coordinates introduces off-diagonal terms in the metric, which for 
the Schwarzschild case (|l8|], pl81) take the form of h{R)X -dX where X are 
the new harmonic Cartesian coordinates. In this work, we prefer to remain 
in isotropic coordinates, where the metric has a simple diagonal form and 
its spatial portion is Euclidean conformally flat, and remove the cross term 
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— V£ • VV> arising from g^T^ by the transformation of the wavefunction, 
Eq. ©. 



A direct substitution of a plane wave- like solution of the form of Eq. (25) 
does not yield the mass shell constraint Eq. (||). This is due in part to 
the existence of the term rj(r) and the fact that the spatial Laplacian V 2 
generates imaginary terms linear in V • p. In the next section we develop a 
WKB solution to Eq. (||). 

5.1 The WKB expansion 



Let us look for a stationary solution to Eq. (31) of the form 



4> (r, t) = Aexp —S(r,t) 

with A constant and where we take the phase to be of the form 

S(r,t) =S(r)-Et. 
Substitution of Eq. (||) and Eq. (||) into Eq. (|l]) yields 

p 2 - (VS) 2 ) + ihV 2 S - h 2 rj(r) = 0, 



(33) 



(34) 



(35) 



where we have used Eq. (||). 

We now expand S(r) in a power series expansion in h via 

oo 
n=0 

Substituting this into Eq. ( j35"l ) yields the set of equations 





0(h°): (VS ) 2 -p 2 
0{h l ) : 2VS ■ VS! 
0(h 2 ) : 2VS ■ VS 2 



iV 2 S 



iV 2 S! - (VS 1 ) 2 + r ] (r) 



n-l 



0(K 



n>?,\ 



2VS ■ VS n 



(36) 



(37a) 
(37b) 
(37c) 



iV 2 S n -! - VS 3 ■ VSn-j (37d) 
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5.2 The eikonal equation 



Eq. (|37a| ) is the dominant contribution to the phase of the wavefunction and 
represents the eikonal equation. We can solve this equation for So via 

So(r) = J r p(x)dx (38) 

where 



E . . E l {n(r)m4\ E n 2 (r) v(r) . 
p=\p\ = —N(r) = -n(rWl- i = ,(39) 



c c V V E J c o c o 

, , E n 2 v(r) , , 

p(r = — • (40) 

co c 

In Eq. ( |39"| ) we have introduced several definitions. First, the magnitude 
p of the spatial momentum is obtained by rearranging Eq. (||). Since the 
energy £ is a constant of the motion for a static metric, the coordinate 
velocity v(r) = dr/dt is a function of position. Its magnitude is given by || 
v(r) = co/n(r) ^Jl — (f2(r)mco/ E) 2 and is related to p as given above. Note 
that for light m = and v(r) = co/n(r) which allows the interpretation of 
n(r) as the index of refraction for massless particles. For massive particles E 
is an extra degree of freedom which specifies the initial speed of the particle. 
The quantity 



N(r)=n(r)\ll-[ nir)mc2 A = n ^< r \ (41) 



E J c 

may be interpreted as the index of refraction for massive de Broglie waves. 
We will return to this point in Section ||. 

So(r) = constant are just the wave front surfaces with normal given by 

VW-rfr, or ^^U_I_, (42) 

where A(r) is the de Broglie wavelength of the particle. Following Holmes 
p2| , let us characterize the wavefronts So(r) by coordinates r = x(l,a,{3) 
where I is arc length along the ray trajectory normal to surfaces of constant 
So, and a and (3 are coordinates used to parameterize the wavefront surfaces 
(for e.g. spherical coordinates). We write the unit tangent vector to the ray 
dx/dl ELS cL vector in the direction of V5o(r) via 

<m = VS (r) = p(r) 

dl \VS (r)\ p ' { > 
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Note we can write 



dSo dx 



dl dl 

Thus, another way to write Eq. ( |38[ ) is 

S (l,a,f3) = 
5.3 The transport equation 



VS = |p(r)| =p. 



pdl 



(44) 



(45) 



We now need to solve the transport equation Eq. ( |37bD , 2VSo • VSi 
iV 2 So- We note for any function F, we can write 



dF/dl = dx/dl ■ VF = 1/pVSo ■ VF 



(46) 



from the ray equation Eq. (|43|). Thus, substituting V5"o • V5i = pdS\/dl 
into the transport equation gives us a first order differential equation for 

Si(l,a,/3), 



dSi _ i 
~df ~ 2p 



:-V 2 S n 



with solution 



Si(l,a,l3) = ~ f dl^. 

2 J p 

A short calculation shown in Appendix [B|, reveals 

V 2 5 d(pj)/dl 



where 



P 



J 



p J 



dx 



d(l, a, (3) 



(47) 



(48) 



(49) 



is the Jacobian of the transformation from the curvilinear ray coordinates 
(l,a,@) to Cartesian coordinates. To prove Eq. ( |49| ) one needs to show (see 
Appendix Irl) 



d,J = JV 



V p 



j dx 



(50) 
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where dx/dl is the unit tangent to the particle's trajectory, normal to sur- 
faces of constant Sq. The net result is that upon substitution of Eq. (|48|) 
into Eq. (f47|), one can perform the integral to obtain 

5i(r(t,a,fl) = Un ( = Ur) (51) 

2 \p{ro)J{r )J 2 

where we define ji{r) = In f ^"jj^j l > an< ^ r o = r |i=o- To the lowest order 

correction O(h ) in the phase and amplitude, we have found the WKB 
approximate solution 



/ $(r) l p(r )J(r ) 
A Vn(r)V p(r)J(r) 6XP 



J p{x)-dx-Et S j 



(52) 



where the subscript 1 on ipi (r, t) indicates that we have carried out the 
WKB expansion to S n= i. Note that Eq. (|5^) does take into account, to 
lowest order, the term (,(r) = ln(fi(r)$ _1 (r) which arises from the covari- 
ant derivative term g^T^. What has been left out are higher order terms 
involving V 2 ^ (r), and V£ (r) in r\ (r), which come from the quantum cor- 
rections terms, S n >2 (which we deal with next). Note that, already this 
solution is valid for strong gravitational fields (as opposed to only weak 
fields as considered by Donoghue and Holstein J5J) with the restriction that 
p 2 S> h 2 r](r), corresponding to 

A(r) {|V£(r)|, |yv^)|} « 1. (53) 

Recall that £(r) = ln(f2(r)<]? -1 (r)) so that the above condition can still 
possibly hold for reasonable distances close to horizon of a black hole, say. 
We investigate this in the next section. 



5.4 Estimation of terms in Schwarzschild metric 

If we had defined <f> (r, t) = u (r) exjp(—iEt/h) in our wave equation Eq. ( |3l| ) 
the result would have been an equation of the Helmholz form 

V 2 u{r) + k 2 (r) (l-^L\ u (r) = (54) 
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where p(r) = Hk(r), k(r) = koN, with ko = loq/cq = E/Ticq. Therefore we 
want to consider the order of magnitude of the terms | V£(r)| 2 /£; 2 (r) and 
\V 2 £(r)\ 2 /k 2 (r) with f(r) = ln(fi(r)$" 1 (r)). 

For the Schwarzschild metric in isotropic coordinates we have 

= 1 ~ 1 'f ) , $ = ^ £ = ln(l-4) J 

1 + 1/P (1 + 1/p) 2 P 2 

n = ^lt 1 //^) ' fc = " §2 ( 55 ) 

where E' 2 = E/(m(%) and p = r/r s , with r s = GM/c 2 , = 0.74 M/M km 
(the the gravitational radius). 

For massless particles, E' — > oo, so that k = k^n and we obtain 

ma^ ■ |Ve(r)|2 - 4 1 1 |V 2 ^)| _ 2 1 {1 + 1/P 2 ) 

k 2 (r) {korsfp^l + l/pf fc 2 (r) ( V,) 2 / (1 + 1/p) 8 

(56) 

Note the above ratios are finite for all values 1 < p < oo (the valid 
range of the isotropic scaled radius p), even though the numerators and 
denominators of the left hand sides of Eq. (156]) each separately diverge as 



p — > 1. From Eq. (|55|), k(p) p _?° k^, so that /co is the usual wavenumber 
of the particle at spatial infinity. For ordinary wavelengths, kor s 3> 1, since 
r s ~ km, so the above ratios remain incredibly small even down to p = 1. 
Even if we were to consider ultra long wavelengths such that kor s ~ 1, the 
ratios in Eq. ( |56"| ) could still be made small for values of p ~ 2, i.e. r = 
2r s , the Schwarzschild radius. In this case, the term rj(r) arising from the 
covariant derivative term g^T^ is essentially negligible compared to k 2 (r) 
all the down to the Schwarzschild radius, and for all intents and purposes, 
our wave equation for massless particles is of the form 

V 2 u(r) + k 2 {r)u{r) = 0. (57) 

However, in the next section it causes no great difficulty to carry terms 
involving rj(r) along formally. 

5.5 Quantum corrections 

The next order h corrections to the phase and amplitude of the scalar wave- 
function arise from Eq. (|37c|). Using Eq. (E(l) we can write a first order 
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equation for dS2/dl whose solution is given by 

fl dl_ 
2p 



S 2 (l,a,(3) 



-S 2 (l,a,(3). 



(58) 



Since £2 is purely real, this is an 0(H 2 ) correction to the phase. 

The remaining equations Eq. (37d) can be formally solved to give 



S n (l,a,(3) = J ^ 



n-l 



(59) 



where the terms in the [ ] , involve only the previously determined quantities 
{S\,S 2 ,... , S n -\} . The first correction to the amplitude is 0(fi 2 ) and is 
given by S3 via 

/' dl 
— [V 2 S 2 - V/x(r) • V5 2 ] = iS 3 (l,a,P). (60) 

Putting this all together, we have to 0(h 2 ) in the phase and amplitude 
r$(r) 

^(r,t) = J_y^3( r ,t) = 



*(*•) p(r )J(r ) 



cxp 



Q(r) y p(r)J(r) 

r 

p(x) ■ dx- Et- h 2 S 2 (r 



exp [-h 2 S 3 (r)] 

, (61) 



where the subscript 3 on ^3 (r, t) indicates that we have carried out the 
WKB expansion to S n= 3. 

Eq. ([3l]) reveals that for spin particles the quantum phase 



s/n = \(fp 



(x)-dx - Et) , 



(62) 



proposed by Stodolsky is only the lowest order (in K) approximation to the 
full phase. This is in stark contrast to the phase of the Dirac wavefunction, 
for which the phase Eq. ( |62|) is exact. The failure of Eq. ( |62| ) to be the exact 
phase for spin particles is directly attributable to the presence of second 
order spatial derivatives in the GRKGE. Since by Eq. (pi]), each component 
of the massless spin 1 field sastifies the GRKGE these remarks also hold 
for the electromagnetic field 4„(i). For the spin values considered in this 
work, only spin 1/2 particles satisfying the GRDE, which contain first order 
spatial derivatives, have the phase of wavefunction given exactly by Eq. 
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6 The Optical-Mechanical Analogy 



In this section we elucidate the optical-mechanical analogy for which the 
path of a particle in a gravitational field can be considered to arise from 



a spatially varying effective index of refraction [H, 1C , [TT|] . To lowest order 
in h for spin 0,1 particles, and exactly for the case of spin 1/2 particles, 
this is the classical path the quantum particle follows. For static metrics in 
isotropic coordinates the magnitude of the momentum p of the particle is 
function of position via Eq. (|3S|), p(r) = (E/cq) N(r). Using this and the 
quantization conditions Eq. @ in the mass shell condition Eq. (||) we can 
define the phase velocity v p h ase and group velocity v g as 

v phase = k=J^Ty (63a) 

V 9 = % =V ^ ( 63b ) 

where the velocity v(r) is given by |J 

v ( r ) = ~T~\ massless particle, (64a) 




— — - J massive particle. (64b) 

Eq. (|63a| ) allows us to identify N(r) defined in Eq. ( |4ll) as the index of 
refraction for massive de Broglie waves. Similarly, Eq. ( |64a ) and the limit 



N(r) n(r) allows us to identify n{r) as the index of refraction for 

massless de Broglie waves. 

We can also add weight to the assertion that N is an index of refraction 
by deriving a geometrical optics ray equation for N. We begin with our ray 
equation Eq. (|43|), multiply through by p and differentiate both sides with 
respect to the arclength I to obtain 

d ( dx \ d dS 

where we have used Eq. (|4|) in the last equality [33]. Substituting in 
(E/cq) N for p on both sides of the above equation yields the geometrical 
optics ray equation for the index of refraction N, 
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Eq. (65) is also directly derivable from the variational principle 34]. 



Ndl = 0. 



(66) 



Further discussion on the interpretation of N(r) as the index of refraction 
for massive de Broglie waves can be found in [p"T|. 



7 The Quantum Phase for Neutrino Oscillations 

As an application, we will calculate the neutrino oscillation formula based on 
the QMP expressions calculated above, for the assumed mixing of massive 
neutrinos following Fornengo et. al. Q. Our interests are two-fold: (1) 
an example illustrating the explicit computation of the QMP and (2) an 
interpretation of the QMP in terms of an effective index of refraction. 

In flat spacetime neutrinos (spin 1/2) produced by the weak interaction 
process are created in a flavor eigenstate \v a ) which is a superposition of 
mass eigenstates i.e. \v a ) = ^2 k U* k Here U is the unitary matrix 
which mixes the different neutrino mass fields. What actually propagates 
is the mass eigenstates, whose energy E^ and momentum p^ are related by 
the mass shell condition E? = (pkCo) 2 + (m^c 2 ,) 2 , and are determined at 
the production spacetime point A. In general, Ej., pk and are different 
for the different mass states. In flat spacetime each of the mass eigenstates 
propagates as \vk(r,t)} = exp(iSfc/ft)|z/fc) where Sk = Pk • x — E^t. 

Neutrino oscillations occur because the different mass states propagate 
differently due to the differences in their energies and momenta. When they 
arrive at a detector located at a spacetime point B which detects flavor 
eigenstates via the weak interaction process, they have developed a relative 
phase shift. Interference between the different mass eigenstates at B pro- 
duces the neutrino oscillations. One assumes the mass eigenstates are pro- 
duced by some coherent process at the spacetime point A and that they are 
detected at the same spacetime point B. The probability that the neutrino 
\v e ) produced at A is detected as at B is given by (for two generations) 
V{y e — > Ufj) = \ {v e \v fJL {B))\ 2 = sin 2 6 sin 2 (5i2/2?i) where 9 is a mixing angle, 
Sk are the phases acquired by the mass eigenstates, and Si 2 = Si — S2. 

Fornengo et al. use the Stodolsky expression for the QMP as given by 
Eq. (ffl), reasoning that this form of the phase is valid independent of the 
particle's spin. From our result Eq. ( |27| ) for spin 1/2 particles, we see that 
this is indeed the correct choice, though not from their original premise. 
Since the neutrinos all begin at the spacetime point A and are detected at 



22 



the spacetime point B, they all experience the same amplitude change as 
given in Eq. (p7|) , so we can ignore it. We write the phase (without h) as 



S 



B 



dt 
dr 



P, 



dr 



(67) 



Fornengo et al. use the following procedure. They compute the phase S for 
a radial light-like trajectory with the mass shell condition given by Eq. (Q). 
This assumes that one is considering ultra-relativistic neutrinos with E\. 3> 
mfcCg. For radial null geodesies we obtain the condition Eq. ( |7a| ) which 
gives dt/dr = n(r)/cQ. The radial momentum p r is given by Eq. (|3^) , 
Pr = (E k /c )N(r) = (E k /c ) [1 - (n^mkcl/Ek) 2 } 1 / 2 . Noting that E k /c 
is a constant momentum along the geodesic we have 



rr B ,t B 



S k = — 




(68a) 



dr. (68b) 



Under the condition of ultra-relativistic neutrinos, the second term under 
the radical is assumed small and the square root can be expanded to first 
order. For the case of the Schwarzschild metric we define p = r/r s and 
carrying out the integral yields 



oSchw 



c ) 2 V E k ) 



\pb-Pa)\ + 



1 

PB 



1 

PA 



+ 



(69) 



This is essentially the form that Fornengo et al. write down in their paper 
(0, Eq. ( |39| ) ) except that here we use isotropic coordinates, and they make 
the further ultra-relativistic approximation Ef. ~ Eq + 0(m\cQ/2Eo) where 
Eq is the energy at spatial infinity for a massless particle. Oscillations then 
occur at phase shifts proportional to (Am^j /2Eq) \ps — pa\ +C(/° _1 ) 
where Amy = rr? k — mj. 

We note two points. First, the gravitational effects are implicit in 
Eq. (|69| ) since p is the scaled coordinate distance. In the presence of gravity, 
the neutrino propagates over the proper distance L p given by L p = J ^Jg rr dr 
= r s J (1 + 1/ p) 2 dp = r s [p — 1/p + 2 Inp]. Second, Eq. ( |68aD allows us to in- 
terpret the phase Sk as calculated in the procedure of Fornengo et al. as the 
integrated "optical path difference" resulting from the difference between the 
index of refraction for a massless particle n(r) and a massive particle N(r) 
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of momentum Ek/co- Note that if the neutrino was massless, N(r) — > n(r) 
and Sk = 0. The authors chose their method of calculation over that of cal- 
culating the phase along the classical trajectory as in Ref. |35| even though 
the final results agree. In the later case, the classical trajectories of different 
massive neutrinos reaching the detection point at the same time must start 
at the production point at different times. Thus, there are initial phases 
for the wave functions that must be added in "by hand." Fornengo's et 
al. approach calculates the interference between mass eigenstates produced 
at the same spacetime point A and detected at the same spacetime point 
B connected by a null geodesies. We see that their "mixed" approach can 
be interpreted as the accumulation of phase due to the difference between 
massless and massive de Broglie waves. 

If the calculations above were repeated for either a massless spin par- 
ticle or say a photon, the derivation would proceed the same except the 
mass would be set to zero in Eq. (|68b| ) . This would imply the lowest or- 
der contribution to the phase, i.e. the "classical phase" Eq. (67), would be 
identically zero. The next contribution to the phase would come from the 
quantum correction terms exp(— iftS^^r)) in Eq. (61) where 52 (r) is defined 
by Eq. pi). 



8 Conclusion 

In this paper we have examined the proposal that the phase of the wave 
function for quantum mechanical particle in curved spactime takes the form 
of Eq. ([[]), as put forth by Stodolsky ||. We investigated the wave equa- 
tions for spin 0,1/2,1 particles in the background of an arbitrary static 
gravitational field which can be written in isotropic coordinates and devel- 
oped explicit plane wave-like solutions. We found that only for the case of 
spin 1/2 does the phase take the form of Eq. (|l]) exactly. This was directly 
attributable to the first order spatial derivative structure of the Dirac wave 
equation. For spin and spin 1 particles the phase takes the form of Eq. (|l|) 
only to lowest order in h, due to the second order spatial derivative structure 
of the corresponding wave equations. We developed a WKB solution for spin 
particles which is also applicable for spin 1 particles. We noted that in a 
gravitational field, the wave function for the generally covariant extension 
of the Dirac equation is not necessarily in addition a solution to the curved 
spacetime Klein-Gordon equation, as is the case in flat spacetime. We find it 
very intriguing that in the presence of a gravitational field, the Dirac equa- 
tion continues to admit exactly a generally covariant extension of a plane 
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wave-like solution with the phase given by Eq. (|l|), while the Klein-Gordon 
and massless spin 1 wave equations only do so to lowest order in fi. 

For the case of spin 1/2 particles we calculated the quantum mechan- 
ical phase appropriate for neutrino flavor oscillations for radial geodesies. 
For spin and spin 1 particles the phase of the quantum wave function is 
predominantly the classical phase as given by Eq. (|I]), with higher order 
quantum corrections. For most applications, especially those involving the 
solar system, the form of the quantum mechanical wave function in curved 
spacetime assumed by Stodolsky would be essentially correct (except that 
the amplitude would vary spatially) for all practical calculations. However, 
for the case of massless spin and spin 1 particles and radial geodesies, we 
showed that the classical phase is zero and the higher order fi WKB phases 
would be the dominant contribution. 

For all the wave equations discussed in this paper we drew an analogy 
for the geodesic path followed by a quantum particle in a static gravitational 
field to motion in a medium with a spatially varying effective index of refrac- 
tion. By examining the momentum of the quantum particle, we were able to 
define an effective index of refraction n(r) and N(r) for massless for massive 
de Broglie waves, respectively. This allows us to extend the classical optical- 
mechanical analogy to the quantum regime for arbitrary static, background 
gravitational fields which can be written in isotropic coordinates. 



A Dirac Gamma Matrices in Flat Spacetime 

In this section we adhere to the notation of the main body of the text and 
use Latin indices {a, b, c} to indicate flat spacetime indices in the range 
{0, 1, 2, 3}. Latin indices in the middle of the alphabet k} refer to spa- 
tial indices {1, 2, 3}. The Lorentz metric is given by r\ a \, =diagonal{l, — 1, — 1, - 
The defining relation for the 4x4 Dirac gamma matrices is 

{7 a ,7 6 }=7 a 7 6 + 7 6 7 a = 2 ?7 ab . (A.1) 

7° is unitary and Hermetian (7 ) 2 = 14x4 = 7° 7° , while the spatial gamma 
matrices 7* are unitary (7*) _1 = (7*)^ and anti-Hermetian (7*)' = —7*, 
with (7*) 2 = — 14x4 = 7*7 1 ' • In this paper we use the explicit representation 

10,11 
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where and I are the 2x2 zero and identity matrix respectively, 7 
{7 1 ,7 2 ,7 3 }, and a = {a 1 , a 2 , a 3 } are the 2x2 Pauli matrices 



a 



1 

1 



a 



-i 

1 



a 



1 
-1 



(A.3) 



The Pauli matrices have the property a 1 cr J = 5^ + i e ljfc a k . Here e ljfc is 
the Levi-Civita symbol with e 123 = 1 and anti-symmetric in all its indices. 
Another useful relationship is (<r • a) (er • b) = a ■ b + icr ■ (a x b) for any 
arbitrary pair of spatial vectors a and b. 

In the derivation of the wave equation for arbitrary spin in Section |3| we 
introduced the anti-symmetric spin matrices 



CT =;l7 .7j = o(T 7-7 7)' 



2 2 
In the representation of Eq. (A. 2) we have 



a ok = i 



a k 



a 



k 



k 

cr fc 



(A.4) 



(A.5) 



B Proof of Eq. ([ 

In this appendix we will prove Eq. 



V 2 S dln(pj) 



P 



dl 



(B.l) 



where J is the Jacobian of the transformation from curvilinear coordinates 
(I, a, (3) (which describe the wave fronts Sq{t) ) to Cartesian coordinates. In 
order prove this relation, we must first prove the following lemma, Eq. (|5C 



B.l Lemma: Proof of Eq. (|50|) 

We want to prove 

'VSo 



d,J = JV 



P 



T _ dx 

JV 'dF 



(B.2) 



which states that the logarithmic derivative of the Jacobian J along a con- 
gruence of ray trajectories is equal to the divergence of the tangent vector 
field of the congruence. 
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A standard result proved in most relativity books (see for e.g. [12] p242- 



243, g7fl p93-94), is that for any matrix ciij, with determinant a, and inverse, 



a 13 = A 3t /a, where A %3 is the signed cofactor of aj,- ; we have 



da 
dx k 



a a 



3 % da ij 



a a dan 



dx k 



aa- 



Let 



J 



dx k 



dx 



for ciij symmetric. 



(B.3) 



(B.4) 



d(l, a, (3) 

be the determinant of the transformation matrix from Cartesian coordinates 
x, to the curvilinear ray coordinates (I, a, (3). Writing out Eq. (|B,3j) with 



etjj — * Jij = dx 1 /dx' 3 , with x % as Cartesian coordinates and x n as curvilinear 
coordinates yields 



dJ 

dx' k 



J 



dx' 3 d f dx 



dx 1 dx' k \ dx' 3 



J 



dx' 3 d 



dx 1 dx' 3 \ dx k 



dx 1 



J 



d ( dx 1 



dx 1 \ dx' k 



(B.5) 

where in the second equality we have interchanged the order of the differen- 
tiations d/dx' k and d/dx' 3 , and have used the chain rule in the last equality. 
If we set k = 1 with x 1 = I, the term in the last () above is just dx/dl, in 
component form. Thus, with k = 1, Eq. ( B.5| ) is just Eq. flB,2| ) in component 
form. 



B.2 Proof of Eq. fl48j): 



Using Eq. (B.2), we want to show 

V 2 S /p = dln{pJ) /dl. 
Using the ray equation Eq. (flgj) in the form 

dx VS (r) 



(B.6) 



(B.7) 



dl p 

we can write the first equality in Eq. (|B.2j ) as (expanding out the divergence) 

din J VSn Vp , V 2 S* _dx , V 2 S dlnp V 2 5 n 

" ~ ~dl 



dl 



+ 



V 



V 



p 



Vlnp + 



P 



dl 



+ 



P 
(Bi 



Solving for V 2 So/p in Eq. (B.8|) yields the desired result Eq. 
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